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ABSTRACT 

Ultraluminous X-ray sources (ULXs) are some of the most enigmatic X-ray bright sources 
known to date. It is generally accepted that they cannot host black holes as large as those 
associated with active galaxies, but they appear to be significantly more luminous than their 
better understood Galactic X-ray binary (XRB) cousins, while displaying an intriguing com- 
bination of differences and similarities with them. Through studying large, representative 
samples of these sources we may hope to enhance our understanding of them. To this end, 
we derive a large catalogue of 650 X-ray detections of 470 ULX candidates, located in 238 
nearby galaxies, by cross correlating the 2XMM Serendipitous Survey with the Third Ref- 
erence Catalogue of Bright Galaxies. The presented dedicated catalogue offers a significant 
improvement over those previously published both in terms of number and the contribution 
of background contaminants, e.g. distant quasars, which we estimate to be at most 24 per 
cent, but more likely ~17 per cent. To undertake population studies, we define a 'complete' 
sub-sample of sources compiled from observations of galaxies with sensitivity limits below 
10 39 erg s . The luminosity function of this sample is consistent with a simple power-law 
of form N(> L x ) oc hyr ' . Although we do not find any statistical requirement for 

a cut-off luminosity of L c ~ 10 40 erg s _1 , as has been reported previously, we are not able 
to rule out its presence. Also, we find that the number of ULXs per unit galaxy mass, S u 
decreases with increasing galaxy mass for ULXs associated with spiral galaxies, and is well 
modelled with a power-law of form S u oc M~ om±om . This is in broad agreement with 
previous results, and is likely to be a consequence of the decrease in specific star formation 
and increase in metallicity with increasing spiral galaxy mass. S u is consistent with being 
constant with galaxy mass for sources associated with elliptical galaxies, implying this older 
ULX population traces stellar mass rather than star formation. 
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1 INTRODUCTION 

The first X-ray imaging observations of galaxies beyond the local 
group, obtained with the Einstein observatory in the early 1980s, 
lead to the discovery of a new population of extremely bright, 
extra-nuclear X- ray sources with luminosities Lx > 10 39 erg s _1 
dFabbianoll 989). Such sources are now commonly referred to as ul- 
traluminous X-ray sources (ULXs), and remain the subject of con- 
tinuing speculation as to how they are able to radiate so brightly. 
The observed luminosities imply some kind of accreting black hole, 
but they are (or at least appear) brighter than the traditional Edding- 
ton limit for the 'stellar mass' black holes in X-ray binaries (BHBs) 
observed in our own galaxy (A/bh~ 10 Mq). However, it is gen- 
erally accepted that these are not extra-nuclear examples of the su- 
permassive black holes (SMBH; A/ B h > 10 5-6 Mq) expected to 
power active galactic nuclei, as such objects should sink to their 
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galaxy c entres comfortably with in a Hubble time due to dynamical 



galaxy centres comfortably withi 
friction llMiller & Colberd2 004). 



Improvements in the available X-ray instrumentation, espe- 
cially with the launch of the XMM-Newton and Chandra obser- 
vatories, combined with multi-wavelength studies have improved 
our understanding of ULXs considerably in recent years, but the 
exact cause of the apparently extreme luminosities has yet to be 
determined unambiguously for any individual source, let alone for 
the population as a whole. Most of the proposed explanations re- 
volve around one of the three following ideas. The first is that 
the black holes powering ULXs are simply much more m assive 
than those in BHBs (see e.g. IColbert & Mushotzkvl Il999h , per- 
haps the long sought after intermediate-mass black holes (IMBHs; 



10 M < M B 



1O 4 M0), accreting at rates substantially be- 



low the Eddington limit as per Galactic BHBs. The second is that 
the black holes hosted by ULXs are essentially the same as those 
observed in stellar mass BHBs, but in an accretion state in which 
they are able to radiate at or even above the Eddington limit (see 



© 2010 RAS 



2 D. J. Walton, T. P. Roberts, S. Mateos & V. Heard 



iRemillard & McClintoclJl2006l for a recent review of the standard 
BHB accretion states). A number of methods by w hich this may be 
physi c ally possible have been suggested (see e.g. JPoutanen et al.l 
120071; iFinke & Bottchej I200I lAbramowicz et alj|l980b . Indeed, 
observational evidence from ULXs now appears to show differ- 
ences between ULXs and Galactic BHBs, suggestiv e of ULXs rep- 
resenting a new supe r-Eddington accretion state dRobertsll2007l ; 
iGladstone etail2009l) . 

The third idea is that ULXs do not emit iso tropically, and ar e 
sources that we view at a favourable orientation dKing et aT]|200lh . 
Depending on the extent to which the emission is anisotropic, this 
can potentially lead to artificially inflated estimates for the lumi- 
nosity of a source, which is usually calculated under the assump- 
tion of isotropy, and again there may be no need to invoke larger 
black holes than those present in stellar mass BHBs. In the con- 
text of ULXs, it has been suggested that the required anisotropic 
emission may be attained through e.g. geometric collimation of ra- 
diation due to a thick accretion disc (possibly the 'slim' discs pro- 
posed to allow super-Eddington emission, see lAbramowicz et al.l 
1980), or through collimation into re lativistic jets similar t o those 
observed in other accreting sources dRevnolds et al.lll997h . How- 
ever, a growing number of ULXs are observed to be embed- 
ded within roughly sph e rically symmetric p hotoionised nebulae 
dPakull & Mirionil 120031 iKaaret et alj 120041) . Detailed study of 
these nebulae suggests they are not supern ovae remnants (SNRs) 
but are apparently inflated by the ULX dPakull & Grisel 120081) . 
Strongly anisotropic emission may be ruled out in these cases 
through morphology and photon counting arguments, so it appears 
unlikely ULXs as a class may be considered 'microblazars'. 

We are now in the position where recent scenarios explaining 
the physical processes governing the high luminosities of ULXs 
combine two or more of the ideas presented above. For exam- 
ple, super-Eddington discs are predicted to inflate and produce 
outflowing winds tha t act to collimate the X-ray emission from 
their centra l regions dBegelman et al.l 2006; Po utanen et al.ll2007l : 
iKind 2008). Also, it has been suggested that a proportion of 
the ULX population may consist of larger stellar remnant black 
holes (up to ~ 80 M0), that may plausibly be formed in low- 
metallicity environments, a ccreting at mildly super-Eddington rates 
dZampieri & Robertsll2009l) . 

As previously stated, unambiguously distinguishing between 
these explanations, notably the IMBH and stellar mass interpre- 
tations, has proven extremely difficult on the bas is of current ob - 
servations. For a recent review on this topic, see iRobertsI d2007t) . 
In large part, this is because reliable dynamical mass estimates 
are not yet available for any ULXs. Until such mass estimates 
become reality, an excellent way in which to enhance our under- 
standing of ULXs is t hrough populati o n stud i es of large ULX sam - 
ples, as compiled by IColbert & Ptakl d2002|) . ISwartz et all d2004h . 
Iliu & Bregmanl d2005h and i Liu & Mirabeld2005l) . For example, 
such work has demonstrated that ULXs are strongly linked with 
recent st ar formation and are prevalently locate d in spiral galaxies 
(see e.g. ISwartz et alj|2004t iGrim m et al. 2003 and Gilfano v et al.l 
l2004ahlLiu et alj|2006l : but especially ISwartz et alj|2009l and refer- 
ences therein). This relationship with star formation strongly sug- 
gests that the majority of ULXs are a form of high mass X-ray 
binary (HMXB), where a massive young star forms the necess ary 
large mass reservoir for the ULX (see e.g. Rappaport et al. 2 0051) . It 
is also suggested that most of the ULXs observed in elli ptical galax- 
ies are likely to be low mass X-ray binaries (LMXBs; Colb ert et al.l 
l2004lHumphrev et alj[2003h . Based on the luminosity function of 
a very large sample of HMXBs, and the relation between the star 



for mation rate and the t otal X-ray luminosity of sta r forming galax- 
ies. lGrimm et al.ld2003l) and lGilfanov et al.ld2004bh argue that there 
may be an upper limit of ~100 Mq to the black holes in standard 
HMXBs, which if true has important consequences for ULXs. 

In addition, the compilation of large samples of ULXs can 
identify extreme examples of these sources. This is pertinent as it 
has become apparent that there exists a rare subgroup of the ULX 
population which display even more extreme X-ray luminosities of 
Lx > 10 41 erg s" 1 . These have bee n dubbed Hyper-Luminous X- 
ray Sources (HLXsi lGao et al.l2003l) . and may be the most promis- 
ing candidates in terms of the search for IMBHs. To date the most 
notable of these is ESO 243- 49 HLX-1, observ ed to reach lumi- 
nosities of 1.2 x 10 42 erg s" 1 dFarrell et al.l2009l). but only a hand 



ful of other candidate H LXs are known (e.g. Davis & Mushotzky 
|2004| ; IJonker et alj|2010l) . Clearly finding more such objects is an 
interesting and potentially illuminating objective. 

Here, we embark upon such a study and present a large cata- 
logue of X-ray detections of candidate ULXs, compiled from Data 
Release 1 of the 2XMM Serendipitous survey (2XMM-DR1, here- 
after referred to as 2XMM). The paper is structured as follows: 
section 2 details the observations used and the process of compiling 
the catalogue, section 3 presents and considers some of the science 
that may be derived from it, and finally section 4 summarises our 
conclusions. 



2 CATALOGUE PRODUCTION 

As stated, the major resource utilised for our sou rce selection 
was th e 2XMM-DR1 catalogue. This is described by Wats on et al.l 
(2009). In brief, this is a catalogue of 246897 detections drawn 
from 3491 public XMM-Newton fields imaged by the European 
Photon Imaging Camera (EPIC) detectors, data spanning the first 
7 years of the mission. It contains 191870 unique X-ray sources 
(as some sources have been multiply imaged, given that XMM- 
Newton has studied some fields on several occasions), and covers 
roughly ~ 1% of the whole sky. Such a large catalogue presents 
a hitherto unrivalled opportunity to study relatively rare source 
classes such as ULXs; in fact, at the modal detection flux for a 
2XMM source (~ 2 x 10~ 14 erg cm " 2 s" 1 in the 0.2-12.0 keV 
band; see Fig. 9 of I Watson et ai] |2009l) one can roughly expect to 
detect all isolated ULXs within the catalogue field-of-view out to 
d ~ 20 Mpc. We therefore took the entire 2XMM catalogue as 
the starting point for our ULX candidate selection. This section de- 
scribes the process of condensing 2XMM into a catalogue of ULX 
candidates, and the properties of the derived catalogue. 



2.1 Initial Correlations and Filtering 

The initial selection of ULX candidates was performed by cross- 
correlating the 2XMM catalogu e with the Third Reference Cata- 
logue of Bright Galaxies (RC3: lde Vaucouleurs et al.1 fl991). This 
catalogue contains 23022 individual galaxies, the majority of 
which have either isophotal diameters D25 > 1 arcminute, to- 
tal B-band magnitude Bt < 15.5 or recessional velocity cz < 
15000 km s _1 . Before cross-correlation with 2XMM, galactic cen- 
tre positions were updated and missing recession velocities ob- 
tained via a cross-correlation of RC3 with the online NASA Extra- 
galactic Database (NED) 1 and the VIZIER/SIMBAD databases 2 . 

1 http://nedwww.ipac.caltech.edu/ 

2 http://cds.u-strasbg.fr/ 
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Figure 1. An illustration of the 2XMM data and ULX selection process. We show data for two galaxy fields: (upper panels) NGC 4472, also including NGC 
4464 (to the north-west of NGC 4472, where north is up in all panels) and NGC 4467 (west of NGC 4472, overlapping its 25 ellipse); (lower panels) NGC 
5194, overlapping with NGC 5195 to its north. In all panels the 25 isophotal ellipse is shown, and the galaxy nuclear position is marked with an asterisk, (left 
panels) A combined EPIC data broad-band (0.2-12.0keV) image with source detections overlaid. The images were extracted as pipeline products from the 
XMM-Newton science archive (XSA), and are convolved with a 1-pixel (= 4 arcsec) HWHM Gaussian mask for display purposes. The images are displayed 
on a logarithmic heat scale, between 2.5-100000 count pixel -1 (NGC 4472), and 1-10000 count pixel -1 (NGC 5194). Point source detections are marked 
by cyan open squares, and extended source detections by green circles with radius equivalent to the measured radius (note that this is capped at 80 arcsecs in 
the 2XMM catalogue), (right panels) Digitised sky survey (version 2) red images of the galaxy fields, with point source positions overlaid. Negative images, 
linearly scaled, are shown for display purposes. The catalogue creation process is highlighted by showing field sources as blue open squares, with sources 
within the 25 great circle but not the 25 isophotal ellipse shown by (thicker) red squares. Sources within the ellipses are shown as cyan open squares with black 
inner lining, with those selected at ULX-like luminosities highlighted by an additional diamond symbol. 
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Galaxies that remained without a recession velocity after these 
measures were excluded from our analysis. Distances to galaxies 
with recession velocities less than 1 000 km s _ 1 , for which peculiar 
motions may dominate the Hubble flow, were obtained by cross- 
correlation with the nearby galaxy catalogue of lTullvl dl988l) . For 
the other galaxies, the Hubble flow is considered a good approx- 
imation, and distances were obtained using a Hu bble constant of 
H = 75 km s" 1 Mpc" 1 to remain consistent with Tullv ( 1988). 

Cross-correlation of RC3 with 2XMM was performed using 
the TOPCAT graphical user interface 3 . An initial match was per- 
formed between 2XMM and RC3 via a conical search, within a 
radius equal to half of the RC3 galactic 25 semi-major axis (where 
the 25 ellipse is denned as the elliptical equivalent to the 25th mag- 
nitude isophote for a galaxy). An elliptical filter was then applied to 
the matched source list to select only those sources that lay within 
the 25 elliptical isophote of all the galaxies, based on the RC3 ma- 
jor and minor axis ratios and position angles. However, a number of 
galaxies included in RC3 do not have position angles listed. In these 
cases, we instead performed a circular match within the minor axis. 
At this point we also discarded all sources flagged as extended, as 
ULXs should by definition be point-like (although a more detailed 
discussion of this exclusion is given in section [231 >. The success 
of this process was verified by visual inspection for a number of 
galaxies; we show the examples of NGC 4472 and NGC 5194/5 in 
Fig.[Tjfor illustration, as well as to highlight the complexities and 
limitations of the 2XMM source selection process. 

X-ray luminosities were calculated for all remaining sources 
using the full band EPIC flux (0.2-12.0 keV), as given in the 
2XMM catalogue, and the distances to the putative host galaxies 
obtained earlier. The 2XMM fluxes are calculated separately for 
each of the three EPIC detectors from the observed count rates, as- 
suming an absorbed powerlaw model for the average source spec- 
trum (T = 1.7, iV H = 3 x 10 20 atom cm -2 ). The quoted EPIC flux is 
then the weighted average of the fluxes obtained for the individual 
detectors dWatson et al.ll2009f) . Although luminosity uncertainties 
are likely to be dominated by the uncertainty in the distance to the 
galaxies, particularly where local peculiar motions may be impor- 
tant, the uncertainties in these distances are not well-quantified (al- 
though any error based on our choice of Ho will at least be system- 
atically propagated through the sample). We therefore calculated 
the uncertainty on the luminosity solely from the uncertainty on 
the flux. As dictated by the standard definition of a ULX, sources 
with Lx < 10 39 erg s _1 were discarded, although those that were 
consistent within 1 a of such luminosities were retained in order to 
include the largest number of possible ULX candidates. We also 
adopted a slightly conservative stance on the quality of detection 
to retain, requiring them to be at least 3.5a detections according 
to 2XMM. This analysis resulted in a large preliminary catalogue 
of source detections at (or above) ULX luminosities, and a small 
sample of slightly fainter X-ray sources. 

2.2 Filtering for Known Contaminants 

The largest identifiable contaminant of our sample of ULXs were 
the active galactic nuclei (AGNs) within the ULX host galaxies. A 
simple cut to remove classic AGNs at Lx 10 42 erg s _1 would 
however be insufficient for our purposes as it would not remove 
the potential low-luminosity AGNs (LLAGNs) contaminating our 

3 Tool for Operations on Catalogues and Tables; 
http ://w w w. star, bris . ac . uk/~mbt/topcat/ 
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Figure 2. Cumulative distribution plot, showing the separation of the fil- 
tered 2XMM detections from the centres of the host galaxies. We plot this 
only for sources in close proximity to the nucleus (r max < 20 arcsecs), 
and show the sample of classic AGNs (solid line), and the comparison 
ULX/LLAGN sample (dashed line) distributions separately. 

sample. Unfortunately, LLAGNs overlap significantly with ULXs 
in terms of observed luminosity, with LLAGN luminositie s as low 
as ~ 10 38 ergs" 1 dGhosh et alj|2008l ; IZhang et alj|2009l) . hence 
making ULXs and LLAGNs indistinguishable by X-ray luminosity 
alone. The only way to circumvent this problem is to remove all 
possible LLAGN candidates; hence we were required to remove all 
sources in close proximity to the nucleus. 

They key question then is: what exclusion radius around the 
centre of the galaxy will optimise the removal of LLAGN can- 
didates, while retaining the maximum number of ULXs? We ap- 
proached this in an empirical fashion by defining a conservative 
maximum separation of each source from the nuclear position, 
r max as the sum of its calculated separation and the 3cr error in 
its position as defined in 2XMM. We then selected the classic 
AGN contamination of the sample by discarding any sources with 
L x < 10 42 erg s~\ leaving 102 likely bona fide AGNs. We plot 
the separation of these sources (quantified as the maximum separa- 
tion, r max ) from their host galaxy central position as a cumulative 
distribution in Fig. [2] The same distribution for all sources with 
10 39 < L x < 10 42 erg s -1 , taken from the same working version 
of the catalogue, is also shown in Fig. [2] in order to illustrate how 
any empirical cut would affect this section of the dataset. 

The initial cut was taken at 5", as this removed 95% of the 
AGN sample, but only 25% of the sources with luminosities in 
the ULX range. As a significant fraction of the latter sample are 
likely to be LLAGNs, the number of actual ULXs discarded is un- 
likely to have been too large. However, an inspection of the remain- 
ing catalogue revealed that there were still a significant number of 
known and probable AGNs present. We therefore tried larger exci- 
sion radii, and found that a value of r- max = 7.5" struck the best 
balance between removing known AGNs and leaving a large sam- 
ple of ULX candidates (< 40% of the ULX-luminosity sample are 
removed, many of which we suspect to be LLAGNs at these small 
nuclear separations). 

The remaining sample was further cleaned with the following 
steps. Firstly, all high luminosity sources (Lx > 10 40 erg s _1 ) 
were inspected by checking their host galaxy nuclear classification, 
their separation from the galaxy centre, and the size of their po- 
sition error. We identified further examples of AGNs that had not 
been filtered by the excision radius due either to a poorly-defined 
nuclear position for the galaxy, or to a large position error, with the 
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Table 1. Some brief details of the detections included in, and the host galaxies contributing to, the presented catalogue and some of its key subsets. Notes: both 
the average X-ray flux fx , and the average total counts per source (summed across all three EPIC detectors) C, are given as the median value across all ULX 
candidate detections. All other averages are the arithmetic mean, and are shown with the calculated error on this value. The average X-ray flux and luminosity 
values are for the full 0.2-12.0 keV energy band. 





Full catalogue 


Spiral sample 


Elliptical sample 


Complete sample 


Total detections 


650 


450 


200 


242 


Individual candidate ULXs 


470 


305 


165 


169 


(with multiple detections 


101 


70 


31 


31) 


Host galaxies 


238 


141 


97 


71 


(with multiple candidate ULXs 


113 


79 


34 


43) 






per host galaxy 




(d) (Mpc) 


37.9 ±2.2 


29.9 ±2.1 


49.6 ±4.0 


13.1 ±0.7 


(Mb) 


-20.5 ±0.1 


-20.3 ±0.1 


-20.7 ± 0.1 


-20.2 ±0.1 






per ULX candidate detection 




/x (erg cm -2 s _1 ) 


4.2 x 10~ 14 


4.8 X HT 14 


3.3 x 10~ 14 


1.4 x 10~ 13 


C(ct) 


293 


314 


279 


1160 


( log L x > (erg s _1 ) 


39.48 ±0.02 


39.46 ± 0.02 


39.54 ±0.04 


39.31 ±0.03 


(HR1) 


0.48 ±0.01 


0.52 ±0.02 


0.41 ±0.03 


0.55 ± 0.02 


( HR2 } 


0.23 ±0.01 


0.26 ±0.02 


0.16 ±0.03 


0.30 ±0.02 


(HR3> 


-0.27 ±0.01 


-0.23 ±0.02 


-0.36 ± 0.02 


-0.22 ±0.02 


(HR4) 


-0.46 ±0.01 


-0.46 ±0.02 


-0.45 ± 0.03 


-0.50 ±0.01 



latter resulting from the source being weakly detected and/or well 
off-axis and/or in a region with a high diffuse background com- 
ponent; these were also discarded. There were also a number of 
high luminosity sources close to (within ~ 15" of) the centres of 
elliptical galaxies that were not classified as AGNs. A conserva- 
tive stance was also taken with these objects, many of which were 
poorly parameterised due to the presence of a strong diffuse ISM in 
the galaxies, as we were concerned the 2XMM catalogue may not 
be differentiating between peaks in the strong ISM and true point 
sources, hence these sources were also excluded. The second step 
related to sources that were multiply-detected in 2XMM, where 
only some detections lay within the r max = 7.5 arcsec excision 
radius (which again was mainly due to poorer position constraints 
on some detections). All detections of these sources were excluded 
from the catalogue. 

These measures remove the contamination of the sample by 
AGNs associated with the host galaxies. However, other AGNs may 
still be present in the field of view, predominantly at cosmic dis- 
tances behind the host galaxies, some of which will be known ob- 
jects. We therefore cross-correlated our remaining catalogue with 
the NED and SIMBAD/VIZIER databases to search for known 
QSOs and other potential contaminants. We found a total of 54 
detections of 39 suspected contaminants coincident with our host 
galaxies, including 17 AGNs, 7 foreground stars and 13 possible 
X-ray detections of recent supernovae within the host galaxies. All 
these detections were discarded from our catalogue. Of course, this 
analysis will not correct for any previously unknown contamination 
of the sample with background or foreground objects, the magni- 
tude of which was relatively large in previous ULX samples. We 
return to this problem in i]2.4l 

Finally, there were a number of other minor issues that re- 
quired attention during the production of the cleaned catalogue. 
Close inspection revealed a number of repeated entries of the same 
detection of individual sources. This occurred where the source po- 
sition lay under the 25 ellipse of multiple galaxies (see Fig Q] for 
some examples). In this situation we assumed the source was truly 
associated with the galaxy whose centre it was closest to, and re- 
moved all other entries of the same detection. In addition, we also 



inspected the combined EPIC images of detections with small sep- 
arations from their respective galactic centres to remove detections 
within the pattern of scattered light produced by the mirror supports 
for bright sources, which are likely to be artificial. 

Once we had our catalogue of good ULX candidates, we re- 
turned to 2XMM to re-append any detections of these sources that 
had previously been discarded. This allowed for sources with vari- 
able luminosity, sources located right at the edge of the 25 isophote 
of their respective galaxies with detection positions scattered on ei- 
ther side, etc. It is after this analysis and filtering procedure that we 
present the catalogue. 

2.3 The Catalogue 

The filtering processes left us with a catalogue comprising of 650 
detections of 470 individual candidate ULXs, 101 of which were 
detected on two or more occasions. Of the detections included, 
121 have L x < 10 39 erg s~\ of which 79 are within \a of the 
classic ULX luminosity limit; the rest are detections of sources 
with variable luminosity which have been observed to emit at 
Lx > 10 39 erg s _1 at other epochs. The catalogue contains 369 
data entries per candidate ULX detection, including the full RC3 
entry for the host galaxy, the full 2XMM entry for the X-ray source, 
and additional information we have added in the derivation of 
the catalogue (distance, X-ray luminosity, etc.). For ease of com- 
parison, we highlight probable common en t ries between o ur cat - 
alogue and those of Icolbert & Ptakl J2002l) Iswartz et all d2004h . 
iLiu & Bregmaril d2005l) and lLiu & Mirabelf d2005h . The catalogue 
contains 367 sources that do not appear to be in any of those works, 
that we flag as 'new' ULX candidates. The position errors on the 
detections are on average only 1", and are always less than 3.5". 
Additional summary details of the catalogue are presented in Table 
[T] as well as a detailed example of some catalogue entries in Ap- 
pendix A. The full catalogue itself is available electronically with 
this paper 4 . 



4 Catalogue to be made available online via MNRAS; available in the 
meantime via private communication 
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Figure 3. Fractional distributions for the absolute B band magnitude (magnitude bins) and the distance (20 Mpc bins) for RC3 galaxies that host ULX 
candidates. Arrows indicate the values quoted in Table [T] Elliptical galaxies are in general brighter, and more distant than their spiral galaxy counterparts. 
For comparison, we also show the distributions of the RC3 galaxies covered by 2XMM-DR1 observations that did not contribute any ULX candidates to our 
sample. Note that the gaps in the histograms are for aesthetic purposes only, and that not all RC3 galaxies have B-band magnitude entries. 



In this section, and much of the rest of this paper, we choose 
to treat the ULX candidates hosted by spiral galaxies and those 
hosted by elliptical galaxies separately. There are very clear rea- 
sons for making this separation. Previous work has shown that the 
X-ray source populations of galaxies are dependent on two fac- 
tors: the mass of the galaxy, which relates to an underlying popu- 
lation of older LMXBs, and the amount of on going star formation , 
which relates to the HMXB population (e .g. IColbert et ai]|2004 
iHumphrev et alj|2003l : iLehmer et al .1 12010). In terms of ULXs, it 
is clear that many of the ULXs in spiral galaxies - particularly in 
hosts with a high star formation rate - are directly related to the 
star for mation and are likel y a type of high-mass X-ray binary sys- 
tem (see lSwartz et al.l2 009 and references therein). This large-scale 
star formation is not present in elliptical galaxies, and so a separate, 
physically distinct type of ULX must be present in these systems, as 
argued bv lKind J2002t) . We therefore distinguish elliptical and spi- 
ral candidate ULXs in our analysis, making the distinction based on 
the numerical index for the Hubble type, T, as registered in RC3. 
Candidate ULXs in types SO/a and later (T 1) are taken to com- 
prise the spiral sample (which therefore includes all morphologi- 
cally irregular and peculiar systems), with the remainder of sources 
in hosts with T < 1 constituting the elliptical sample (which in- 
cludes both elliptical and lenticular host galaxies). 

Table Q] shows that we detect a factor of ~2 more ULX 
candidates in spiral galaxies than in ellipticals, and we have a 
slightly higher proportion of multi-epoch datasets for spiral can- 
didate ULXs. We also have a factor ~2 more host spiral galaxies 
than ellipticals, although as shown in Fig.|3]the elliptical hosts have 
a mean distance ~ 70% greater than spirals, and also are on aver- 
age bigger and brighter with a mean absolute magnitude nearly half 
a magnitude higher. The mean log X-ray luminosity of the candi- 



Table 2. Energy ranges covered by the 2XMM energy bands. 



2XMM Band 


Energy 




(keV) 


1 


0.2-0.5 


2 


0.5-1.0 


3 


1.0-2.0 


4 


2.0^1.5 


5 


4.5-12.0 



date ULX detections in each sample are approximately the same 
(see Fig. |4); however, given the greater distance to the ellipticals 
the candidate ULXs possess a fainter median X-ray flux. Despite 
this, we accumulate similar median numbers of X-ray photons per 
detection in both samples (also demonstrated in Fig. |4), since the 
elliptical galaxies typically have longer observations. Note that the 
values quoted in Table Q] are calculated for all the detections in- 
cluded in the catalogue, if we limit ourselves to detections above 
10 39 erg s , the mean log(Ljf) values increase to 39.61 ± 0.03 
and 39.63 ± 0.04 for spiral and elliptical detections respectively, 
which are extremely similar. 

Perhaps the most interesting distinction between the samples 
comes in terms of the 2XMM X-ray hardness ratios, defined as: 



HR 



Cb — Ca 



(1) 



Cb + Ca 

where Ca and Cb are the count rates of the two bands under con- 
sideration. The energies covered by the five 2XMM energy bands 
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Figure 4. Fractional distributions for the X-ray luminosities and the number of photons detected from the ULX candidates (0.5 dex bins) included in the 
presented catalogue. Arrows indicate the values quoted in Table [JJ The mean X-ray luminosity and the median numbers of counts received are similar for 
detections associated with spiral and elliptical galaxies. 
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Figure 5. Fractional distributions for the 2XMM EPIC hardness ratios HR2 and HR4 (AHR = 0.2 bins) for the ULX candidates included in the catalogue. 
Arrows indicate the values quoted in Table[JJ Detections from spiral galaxies appear harder at low energies than their elliptical counterparts, but this discrepancy 
is not seen at higher energies. This may be an indication that the sources in spiral galaxies are more absorbed. 
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are given in Table|2] HR1 combines bands 1 and 2, HR2 combines 
bands 2 and 3, etc. The 2XMM hardness ratios are normalised to be 
in the range [-1,1] with harder spectra having more positive hard- 
ness ratios. Although Fig. [5] shows the detections display a wide 
range of hardness ratios, the mean hardness ratio values are con- 
strained to about ±0.02 in each case; hence the differences between 
the spiral and elliptical sample are significant at the ~ 3 — 4.5a 
level for hardness ratios HR1 , HR2 and HR3. In each case the spiral 
hardness ratio indicates the sources are somewhat spectrally harder 
than their elliptical counterparts, with this spectral distinction not 
apparent in HR4. To demonstrate this we show the number distri- 
butions of HR2 and HR4 in Fig.[5]for both the elliptical and spiral 
galaxy detection populations. The difference in the mean values for 
HR2 is largely driven by an absence of detections with hard spectra 
below 2 keV from elliptical galaxies. We suggest this may primar- 
ily be due to the presence of more neutral gas in spiral galaxies, par- 
ticularly in the star-forming regions hosting the ULXs, which will 
lead to more absorbed and hence harder spectra. However, a higher 
column physically associated with each ULX in spiral galaxies - 
perhaps a result of mass-loss from the likely high-mass stellar sec- 
ondaries in such systems - is also an intriguing possibility. What- 
ever the cause, this difference supports the decision to separate our 
candidate ULXs into two distinct samples. 

Within the catalogue there are also 5 sources that appear to 
display HLX luminosities, none of which have previously been 
included in ULX catalogues to date. The brightest of these is 
2XMM J134404.1-271410, which was observed at Lx = (2.82 ± 
0.24) x 10 41 erg s _1 . However the closest HLX included in our 
catalogue, and hence potentially the easiest to study, is 2XMM 
J01 1942.7+032421 (hereafter NGC470 HLX1), associated with 
the relatively nearby (d ~34 Mpc) spiral galaxy NGC 470. This 
source was observed to reach a luminosity of Lx = (1.53±0.08)x 
10 41 ergs -1 . These sou rces will be the subject of forthcoming 
work dSutton et alfcOld : Sutton et al. in prep.), which will include 
recently obtained follow-up observations of NGC 470 HLX1. 

2.3.1 The complete sub-sample 

In order to undertake certain statistical studies of the population, 
it is necessary to define a 'complete sample', i.e. the subset of the 
catalogue compiled from observations of galaxies in which all the 
ULX candidates present within that galaxy should be detected. By 
doing so, we avoid including an artificial bias for brighter sources. 
Note that in this context we consider an observation of a galaxy to 
be complete if it is possible to detect all the sources that match our 
selection criteria; this does not address the limitations with respect 
to completeness inherent within these criteria, which are discussed 
in S 

To determine whether an observation of a galaxy can be con- 
sidered complete there are two limiting flux sensitivities that must 
be considered and compared. The first is the sensitivity of the ob- 
servation, fobs, i.e. the minimum count rate that a source must have 
to be detected at a certain position in the detector with a given de- 
tection significance. This value is mainly determined by the local 
effective exposure and background level. In order to compute a sen- 
sitivity map for each XMM-Newton observation in the soft (0.5-2.0 
keV) and hard (2.0-12.0 keV ) bands we used an e mpirical approach 
( see ICarrera et al.ll2007l and iMateos et al]|2008l for a detailed de- 
scription of the method). Briefly, if Poisson statistics hold, it is pos- 
sible to determine the minimum count rate that a source must have 
to be detected with a given significance. This Poissonian count rate 
is calculated from the total number of background counts and mean 
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Figure 6. Number distributions of the Hubble type for all galaxies with 
complete observations, and galaxies with complete observations that host 
one or more ULX candidates. It is clear that ULXs are preferentially located 
in early type spiral galaxies. 

effective exposure around the source position. However, source de- 
tection and parameterisation in the 2XMM pipeline is a compli- 
cated process that involves a simultaneous maximum likelihood fit 
of the distribution of source counts detected with each EPIC cam- 
era in different energy bands convolved with the telescope's point 
spread function (PSF). Therefore the count rate values resulting 
from the source detection process are bound to deviate from the 
as sumption of pure Po isso nian statistics. Howe ver, as it is shown 
in lCarrera et al.l (120071) and lMateos et al.l ( 120081) . there exists a lin- 
ear relationship between Poissonian count rates and those derived 
from the XMM-Newton pipeline. This relationship can be used to 
empirically correct the pipeline count rates from all non Poissonian 
effects. This approach is used by the Flux Limits from Images from 
XMM-Newton (FLIX) server provided by the XMM-Newton Survey 
Science Centre (SSC) 5 . Throughout this work, we utilise sensitivity 
maps calculated for a 5a detection significance. 

The second is the flux of a source emitting at 10 39 erg s _1 
at the distance of the galaxy in question, which we refer to as 
the ULX luminosity limited sensitivity, f„i x . In order to consis- 
tently compare fobs and f u lx, we calculate f u i x for the full XMM- 
Newton bandpass (0.2-12.0 keV), then distribute the sensitivity be- 
tween the hard and soft bands according to an average ULX spec- 
tral shape. We adopt a simple absorbed power-law model, with 
(r) = 2.2 and (7V H ) = 2.4 x 10 21 atom cm " 2 , determined from 
the sample presented in Gladstone et al. U200l . Using PIMMS we 
find that approximately 36 per cent of the flux from ULXs should 
be observed in the soft band. 

Having calculated f u i x for each energy band, we compare this 
with fobs for each observation of each galaxy across its projected 
area (excluding chip gaps, etc., as well as the circular region of 
radius 7.5" excluded from our analysis during our consideration 
of nuclear sources). This comparison is performed separately for 
the EPIC-pn and two EPIC-MOS detectors. An observation of a 
galaxy in a particular band with a particular detector is considered 
complete if f u i x is greater than fobs for the whole of the considered 
galaxy area. We take an observation of a galaxy to be complete in 



http ://w w w. ledas . ac .uk/flix/flix .html 

6 http://heasarc.nasa.gov/docs/software/tools/pimms.html 
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general if it can be considered complete in either of the bands for 
any of the detectors, as regardless of whether it is complete for 
only one band in one detector or both bands in all detectors, all 
the ULXs detectable with our selection criteria should appear in 
the catalogue. On investigation, we find that any observation that 
is considered complete in the hard band is always considered com- 
plete in the soft band as well. 

The complete sub-sample consists of 242 detections of 169 
discrete sources, including 130 sources which are observed to radi- 
ate at or in excess of 10 39 erg s" 1 . Of these, 106 are found in 52 
spiral galaxies, and 24 in 12 elliptical galaxies. The spatial coverage 
of 2XMM includes complete observations of a total of 108 spiral 
and 56 elliptical RC3 galaxies (i.e. galaxies for which the com- 
pletion criteria are met, regardless of whether they actually host 
a ULX candidate). In Fig. [6] we show the number distributions of 
the Hubble type for all the galaxies with complete observations, as 
well as those with complete observations that host ULX candidates. 
It is clear that ULXs are preferentially located in early type spiral 
galaxies, as would be expected give n the association be tween ULXs 
and regions of high star formation dSwartz et al.l2009h . The lack of 
ULX candidates in elliptical galaxies is a consequence of their rel- 
atively low star formation rates, while the lack of ULX candidates 
in late type spirals arises due to the majority of these being dwarfs 
galaxies, which do not contain enough stellar mass to host ULXs. 
Of the 169 sources included in the complete sample, 106 are 'new' 
detections, of which 71 are observed to radiate at or in excess of 
10 39 erg s _1 . Finally, there are also 39 sources included that are 
not observed to radiate at Lx ^ 10 39 erg s _1 at any epoch con- 
tributing to the complete sample, although 8 of these sources are 
observed at ULX luminosities at other epochs that contribute only 
to the general catalogue. These fainter sources are also primarily 
observed in spiral galaxies, with only 7 of the 39 associated with 
ellipticals. 

In Table[JJand Figures|4]and|5]we provide a comparison of the 
complete subsample and the catalogue as a whole. As expected, the 
mean distance to host galaxies for the complete sample is signifi- 
cantly smaller than that of the catalogue as a whole (all galaxies 
with complete observations are within 24 Mpc), and the host galax- 
ies are fainter on average. In addition, although the detections in 
the complete sample typically have many more counts, their av- 
erage luminosity is lower than that of the whole catalogue. This 
arises naturally from the fact that a large portion of the detections 
in the full catalogue arise from observations of galaxies that did not 
have the sensitivity to detect sources at 10 39 erg s - 1 . Consequently, 
there is a slight bias towards brighter sources in the full catalogue, 
which is not present in the complete subsample. Finally, there are 
also some slight differences in the average hardness ratios between 
the complete sample and the full catalogue, which will be discussed 
in more detail in forthcoming work (Walton et al. in prep.). 

2.4 Estimating the Contribution from Unknown 
Contaminants 

When producing a catalogue of any specific astronomical object it 
is important to consider what fraction of the sources included are 
likely to be undesirable contamination, e.g. background quasars in 
this case. Here we attempt to quantify the expected contribution 
from such unknown contaminants. In brief, this requires estimation 
of the number of background sources that should be observed at 
some minimum detection significance (for all RC3 galaxies ob- 
served), and comparing this with the actual number of sources 
detected at or above that minimum significance. Like the consid- 



eration of observation completeness, our calculations have been 
performed separately for each EPIC detector in the soft and hard 
bands, as dictated by the availability of the sensitivity maps. 

Since the catalogue compiled is of sources with Lx > 10 39 
erg s~ , only background contaminants with apparent luminosi- 
ties also in this range needed to be considered, so we again make 
use of the full band fuix calculated previously. However, we are 
now considering a different class of source, most likely back- 
ground quasars, which on average may not have their flux dis- 
tributed between the so ft and hard bands in the same way as ULXs. 
IPiconcelli et al.l (2003) find that observed spectra of the majority 
of background quasars are adequately modelled with a powerlaw 
continuum, modified by Galactic absorption, so we again adopt a 
simple absorbed power law model for the average spectral shape of 
background sources. In this calculation we use t he average pho- 
ton in dex of (r) = 1.59 ± 0.02 obtained by IPiconcelli et all 
d2003l) and the Galactic col umn density in the direction of each 
galaxy dKalberla et afl fcoOS). For each galaxy the fraction of the 
0.2-12.0 keV flux that should be observed in the soft and hard 
bands (0.5-2.0 and 2.0-12.0 keV respectively) was determined us- 
ing PIMMS, and f u i x was recalculated for each band accordingly. 
On average ~30 per cent of the flux from these sources should 
be observed in the soft band. We again compare these sensitivities 
to the observation sensitivities provided by the sensitivity maps, 
fobs' The true limiting sensitivity fu m relevant to this calculation 
is whichever of these two is greater, as expressed by equation|2] 

j? J fuix fuix ^ fobs 

yjobs Julx \ Jobs 

Adopting this form simultaneously accounts for observations of 
galaxies deep enough to detect sources fainter than ULXs, and ob- 
servations in which not all the ULXs in a galaxy could be detected. 
Using equation [2] modified sensitivity maps of fu m were gener- 
ated. 

It has long been known that the majority of the Cosmic X- 
ray Background (CXB) can be ac counted for by e mission from 
discre te sources, as summarized in Hasingej] d2004l) . iMoretti et al.l 
d2003l) investigated the number of such sources that should be re- 
solved (N(>S), per unit sky area) as a function of flux sensitivity, 
S, finding that this was well modelled as a smooth broken power 
law in both the hard and soft bands (note that the hard band used in 
IMoretti et alj|2003l only covers 2.0-10.0 keV, however we expect 
that any corrections due to the slight differences in energy range be- 
tween their hard band and the hard band used here will be minimal, 
so it should still be acceptable to compare the two). These relations 
were used to convert our limiting sensitivity maps into maps of the 
number of background contaminants expected from each pixel (one 
pixel covers a 4 x 4 arcsec area of the sky). The total background 
contribution from each galaxy was taken as the sum of all the indi- 
vidual pixel contributions that fell within their 25 isophote (again 
excluding the nuclear region, chip gaps, etc.), and the total back- 
ground contribution expected in the catalogue was then the sum of 
the contributions from all the RC3 galaxies covered by observations 
included in 2XMM. Where galaxies have been observed on multi- 
ple occasions, only the contribution from the deepest observation 
was considered. 

In order to compare the number of expected background 
sources with the number included in the catalogue we generated 
hard and soft band maximum likelihoods (MLs) using those given 
for the five basic energy bands included in 2XMM. We make the 
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Table 3. Fractional background contamination estimates calculated for the hard and soft bands for the PN and MOS detectors (see text). The analysis is also 
broken down into separate considerations of the catalogue subsets from spiral and elliptical type galaxies, as well as the complete sub-sample. 



Estimated Contamination (%) 




Soft Hard 


Soft Hard 


Soft Hard 




All Galaxies 


Spiral Galaxies 


Elliptical Galaxies 


PN 


33.9 ±4.6 12.9 ±3.7 


28.3 ±4.8 10.8 ±3.5 


46 ±10 19 ±10 


MOS1 


iJiJ.O _1_ •J • -L _LC7.i_> _1_ rt . Zj 


oq q 4- z x. i o s + 4 n 


rt 1 _1_ J__L ij C7 _l_ IZi 


MOS2 


37.1 ±5.4 18.5 ±4.1 


28.8 ±5.3 12.5 ±4.0 


58 ±14 35 ±11 


Weighted 
Average 


24.0 ±4.1 


17.9 ±3.8 


39 ± 5 




All Complete Galaxies 


Complete Spiral Galaxies 


Complete Elliptical Galaxies 


PN 


32.9 ±6.6 14.0 ±4.6 


28.0 ±6.5 11.8 ±4.5 


54 ±21 23 ± 15 


MOS1 


37.2 ±7.3 17.9 ±5.0 


31.9 ±7.3 14.1 ±4.9 


58 ±22 34 ± 17 


MOS2 


35.8 ±7.1 17.5 ±5.0 


30.5 ±7.1 13.8 ±5.0 


57 ±22 32 ±16 


Weighted 
Average 


22.5 ± 4.4 


18.5 ±4.0 


39 ± 6 



approximation that the hard and soft ML values are the sum of 
the likelihoods of the basic band they each encompass, which is in- 
creasingly good for higher detection likelihoods. We then counted 
the number of sources (note that we counted sources rather than de- 
tections, as we only considered the deepest observation for galaxies 
observed multiple times) with a detection significance of at least 5a 
(ML > 15). These were compared with the number of estimated 
background sources calculated previously to obtain fractional con- 
tamination estimates for each band and each detector, after taking 
into account the known background contaminants already removed 
(see Q2.2\ . Our estimates are given in Table [3] with quoted uncer- 
tainties based on counting statistics. We also present the results of 
the same calculation for the complete sub-sample of the catalogue, 
as defined in i]2.3.1| 

The average contamination is found to be ~24 p er cent, sim- 
ilar to that estimated for previous ULX catalogues dSwartz et al.l 
120041) . It is clear from Table |3] that the subset of detections as- 
sociated with elliptical galaxies should be significantly more con- 
taminated by unidentified background sources than the detections 
associated with spiral galaxies. Such a result may have been ex- 
pected given the known observational characteristics of the ULX 
population which, as highlighted in i]2.3l are also reflected in this 
catalogue. ULXs a re often associated with star forming regions 
dSwartzetal.ll2009l) , which are more common to spiral galaxies 
than elliptical galaxies, hence ULXs are more commonly observed 
in spiral galaxies. In addition, the observed elliptical galaxies are 
in general both larger and further away than the observed spiral 
galaxies, hence covering a greater cumulative area on the sky at 
lower ULX-detection fluxes. This combination leads naturally to a 
higher estimation of the fractional contamination. 

Finally, we note briefly that the average contamination esti- 
mates presented in Table[3]may be better considered as upper lim- 
its to the fractional contamination rather than true predictions of the 
likely value. This is because our calculation does not correct for the 
flux extinction that will be suffered by background sources due to 
the gas and dust in the galaxies these sources are falsely associated 
with. Such considerations are especially relevant for spiral galax- 
ies, which should contain higher columns of neutral gas and dust 
than their elliptical counterparts, as suggested by the hardness ratio 
comparison discussed in section |2~3l As such, although we quote 
values based on the calculations for both the soft and hard bands to 



be conservative, the contamination estimates for the hard band, in 
which effects due to absorption will be greatly reduced, are likely 
to be more representative of the true contamination. This is there- 
fore likely to be closer to ~ 17 per cent for the whole catalogue, and 
as low as ~12 per cent for the spiral galaxy population. 

2.5 Limitations of the Catalogue 

The largest dedicat ed catalogue of U LX candidates compiled to 
date is presented in lSwartz et al.l I T2004I) . and includes 154 discrete 
sources with an estimated 25 per cent background contamination. 
The catalogue presented here, with 475 sources and at most a simi- 
lar contamination represents a significant improvement in number ' , 
and should provide a useful resource for the ongoing investigations 
into the properties of the ULX population; see Sj3]for examples of 
possible analyses that may be performed. However, before embark- 
ing upon any such study it is important to consider the limitations 
of this catalogue. 

The fractional contamination could be further reduced by 
highlighting and removing additional contaminants with multi- 
wavelength observations and correlations. Such analysis is beyond 
the scope of this paper, but coul d be performed in future work in- 
volving follow-up observations. Won g et all d2008l) present an ex- 
ample of highlighting background AGNs with optical observations, 
but spectra from other wavelengths could also be used to detect 
characteristic features of AGNs, e.g. broad and narrow emission 
lines from Seyfert galaxies. 

Perhaps the major limitation on this catalogue is its incom- 
pleteness. There will inevitably be a population of ULXs residing 
in the galaxies observed that are not included in this catalogue. Ob- 
viously the observations of some galaxies were not deep enough 
to detect all the sources down to 10 39 erg s -1 . In addition, a num- 
ber of legitimate ULXs will have been discarded in the empirical 

7 Since this work was submitted, a large catalogue of gen eral X-ray point 
sources in nearby galaxies has been published by Liu ( 201 1), compiled from 
Chandra observations. That catalogue contains 300 candidate ULXs within 
the D25 galaxy regions, and an additional 179 candidates in the region be- 
tween the D25 and 2D25 ellipses, although we caution that a large fraction 
(~60 per cent) of the latter sample are likely to be contaminats. 
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nuclear cut and by discarding extended sources. In terms of the lat- 
ter reason, star forming regions often appear as diffuse, extended 
emission and due to their well documented association with these 
regions it is likely some ULXs will be embedded within them. 
With its fairly modest spatial resolution, XMM-Newton can have 
difficulty resolving ULXs in this situation, so by excluding ex- 
tended sources we are also excluding any ULXs that may be 'hid- 
den' within them. A particularly good example of this problem for 
XMM-Newton is the relatively nearby (d ~ 35 Mpc), intense star 
forming galaxy NGC 3256. XMM-Newton observations reveal an 
extended X-ray emission regi on centred on the ga laxy, with a soli- 
tary ULX outside this region ferritins et aI1l2004l) . whereas Chan- 
dra has resolved the central regions to reveal a populatio n of ULXs 
embe dded within regions of strong diffuse emission dLira et all 
2002). In addition, detections that appear as extended with XMM- 
Newton may actually be emission from an unresolved population 
of point sources. This may even be the case for a number of detec- 
tions that appear as point sources with XMM-Newton, a possibility 
which becomes increasingly likely for more distant galaxies. The 
combination of these effects is difficult to quantify, but the net re- 
sult is likely to be that there are additional ULXs in some of the 
galaxies observed that are not included in this catalogue, even for 
galaxies that have observations we have considered as 'complete'. 
Observations of these galaxies with Chandra should go some way 
to addressing these issues, owing to its significant improvement on 
the angular resolution of XMM-Newton. 



3 EXAMPLE APPLICATIONS OF THE CATALOGUE 
DATA: ANALYSIS AND DISCUSSION 

In order to demonstrate the scientific potential a resource such as 
this catalogue offers, we now present a brief analysis highlighting 
some characteristics of the derived population based on an analysis 
of the statistically complete sample. We note that the main advan- 
tage enjoyed by the 2XMM catalogue over comparable resources 
is the higher photon statistics per detection, enabling more detailed 
studies of source characteristics over a broader energy range. This 
will be explored in future work (Walton et al. in prep.). 



3.1 Luminosity Function 

One way in which we may investigate the ULX population is 
through construction of its luminosity function. In order to avoid 
the estimation of corrections to account for incompleteness, we 
make use of the previously defined complete subset of the pre- 
sented catalogue, for which such corrections should not be nec- 
essary, or at least negligible (even these data may suffer from some 
level of incompleteness as we are limited by the XMM-Newton an- 
gular resolution; see §2.5\ . There are a number of galaxies con- 
tributing to the complete sample that are found to have multiple 
observations that can be considered complete. Although there is 
some evidence that X-ray luminosity functions don't change sig- 
nificantly from e poch to epoch, ev en if the sources contributing are 
variable (see e.g. Zez as et al . 2007), we adopt a conservative stance 
and only consider the longest observation of each galaxy. 

These observations contribute 121 discrete sources radiating 
at or above 10 39 erg s" 1 . We group these into equally spaced 0.1 
dex logarithmic luminosity bins, with the uncertainty on the num- 
ber of sources in each bin given by Poisson counting statistics. Us- 



ing the SHERPA modelling package , we fit the luminosky func- 
tion by minimizing a form of the Cash statistic (C; lCastJll979l) 
adapted to be an analogue to % 2 (i- e - a reduced statistic of C v ~1 
for acceptable fits), and find that the data are well modelled with a 
single power-law up to the highest luminosity sources in the com- 
plete sub-sample (Lx — 5.5 x 10 40 erg s _1 ); the differential and 
integral forms of such a relation are given in equations [3]and[4] 



dN 
dL x 



oc L\ 



N(> L) oc L~ 



(a-l) 



(3) 



(4) 



This simple model gives a reduced statistic of C v = 0.89, with the 
best fit value of the exponent a = 1.96 ± 0.11 (here and throughout 
this section quoted uncertainties are the la uncertainties on one 
parameter of interest). On investigation, we find that the quality 
of fit, and the exponent obtained are not dependent on our choice 
of luminosity bin size. In Fig. [7] we show the integral form of the 
luminosity function (equation[4), rebinned so that there are approx- 
imately 5 additional sources from one luminosity bin to the next for 
clarity, and the best fit power-law model. The bin centroids in Fig. 
[TJare the average of log(Lx) for the additional sources included in 
each bin. 

We also investigate the luminosity functions for the spiral and 
elliptical galaxy ULX populations separately. Of the 121 sources 
considered previously, 97 are located in spiral galaxies and 24 in 
elliptical galaxies. Again, we apply a simple power-law model to 
each, and find that this provides an excellent representation of the 
data in both cases (C„ = 0.97 and 0.7 respectively), with a sp i ra i 
= 1.85 ± 0.11 and tt e Ui P ticul = 2.5 ± 0.4. As with the full com- 
plete sample shown in Fig. [7] in Fig. [8] we show the integral forms 
of the spiral and elliptical luminosity functions, rebinned in a sim- 
ilar fashion. Although the exponents for spiral and elliptical galax- 
ies agree within their 2o uncertainties, due to the large uncertainty 
on the elliptical exponent, which in turn is a consequence of the 
relative lack of sources from those galaxies, there is a tentative 
suggestion that the luminosity function for sources in elliptical 
galaxies is steeper than for those in spiral galaxies. This would 
be consistent with the steep ening of the elliptica l lu minosity func- 
tion se en in the samples of ISwartz et alj d2004l) and lColbert et al.l 
d2004l) . Such a result would not be surprising given the expecta- 
tion that the sources in elliptical galaxies are likely t o be LM XBs, 
while those in spiral galaxies are likely to be HMXBs dColbert et al.l 
2004; Hum phrey et al.l l2003). and the shapes of the 'universal' lu- 
minosity functions com piled for these two classes of BHB, see 
iGrimm et al. | d2003h and Gilfanov ( 20041) . The exponents obtained 
here are consistent within uncertainties with the relevant exponents 
for the HMXB and LMXB luminosity functions obtained in those 
works, both of which extend to luminosities below 10 39 erg s _1 . 
This might suggest that the same process, or combination of pro- 
cesses, responsible for forming the lower luminosity BHB popula- 
tion are also responsible for forming the majority of ULXs. 

It has previously been reported that there may be a high lumi- 
nosity cut-off in the ULX lum i nosity function at Lx ~ 2 x 10 4 



erg s * , s ee e.g. Swartz et al. J 20041) but especially IGrimm et al.l 
(2003) and lGilfanov et al.ld2004tC who argue strongly for an upper 
limit to the luminosity of HMXBs. Given the large overlap between 
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Figure 7. The luminosity function derived from the complete subset of the 
catalogue. The best fit single power-law model (see text) is also shown. 



the ULX and HMXB populations this is also relevant for the ULX 
luminosity function. Such a cut-off may arise due to there being an 
upper limit to th e mass of the blac k holes HMXBs are expected to 
contain, which Grim m et al estimate to be ~ 100 Mq under 

the assumption that the Eddington limit is not violated. Although 
the data are already reproduced well with a single power-law func- 
tion, we also investigate a broken power-law model, as expressed 
in equation [5] 



dN 
dL x 



T — «1 

L x 
L x 



Lx < L c 
L x > L c 



(5) 



However, this more complex model provides a negligible improve- 
ment in the fit to the data: we find AC = 0.1 for 2 extra degrees of 
freedom, the break luminosity L c is not well constrained (although 
its best fit value is 1.1 x 10 erg s -1 , similar to previously proposed 
values for L c ), and the two expon ents are cons i stent within their 
la uncertain ties. Since the work of lGrimm et al] (f2003) focuses on 
HMXBs and lSwartz et alj d2004 ) only see the cut-off in the lumi- 
nosity function of spiral galaxies, we repeat the comparison using 
only the sources from complete spiral galaxies. Unsurprisingly, as 
our overall luminosity function is dominated by sources from spi- 
ral galaxies, we find exactly the same results, with an improvement 
now of AC = 0.2 for 2 extra degrees of freedom. 

We must remind and caution the reader that our luminosity 
functions are generated from sources from ensembles of galax- 
ies. Although the broad, general shape should remain robust, non- 
systematic errors in the estimated luminosities that differ from 
galaxy to galaxy due to e.g. poorly estimated galaxy distances that 
could smear out complex features such as a cut-off. If present, 
the effect of such uncertainties would be especially prevalent at 
high luminosities where there are fewest sources. In this work we 
have adopted a simple Hubble Flow distance for galaxies with 
cz > 1000 km s" 1 and the distance quoted in iTullvl ( fl988h for 
galaxies with cz ^ 1000 km s~ 1 ; a more careful analysis of galaxy 
distances, which is beyond the scope of this work, might clar- 



ify this issue. Other biases may also be present too, for example 
background contamination. Removing these sources might alter the 
shape of the pure ULX luminosity function. However, a repeat of 
the work in i]2.4l for sources with Lx above 5 x 10 40 erg cm -2 s _1 
(limited to within 100 Mpc) shows ~30 per cent of the sources 
at these luminosities in the catalogue are likely background con- 
tamination, and only ~15 per cent of the spiral galaxy population 
(Sutton et al, in prep.). This is fairly similar to that found for the 
whole ULX luminosity range, so it seems that the fractional con- 
tamination remains roughly constant with luminosity, hence remov- 
ing such sources would primarily effect the normalisation rather 
than the shape of the luminosity function. 

The luminosity functions presented are still ultimately lim- 
ited by low number statistics at the highest luminosities, which 
severely limits the comparison of complex models. So, although 
we conclude that our data do not require anything more than sim- 
ple power-law models, we cannot rule out the presence of an in- 
trinsic cut-off in the ULX luminosity function at L c ~ 2 x 10 40 
erg s~ . However, the complete sub-sample has 5 sources in its 
highest bin, which nominally lies above this break. If the value 
adopted for Ho is significantly underestimated then source lumi- 
nosities might be systematically overestimated, but this seems un- 
likely gi ven the currently favou red value of Ho = 72 ± 8 km s _1 
Mpc -1 iFreedman et alj|200ll) . Given the expected level of con- 
tamination it is also unlikely these are all foreground/background 
sources. They must therefore be regarded as bona fide ULX can- 
didates beyond the previously reported limit, which must call into 
question the veracity of this proposed feature. Indeed, the full cata- 
logue contains an additional 38 sources above this proposed upper 
limit, although we note again that Chandra follow-up of all of these 
high luminosity sources is strictly required to confirm their nature 
as single point sources. 

Even though these additional sources may not be included in 
our quantitative analysis of the luminosity function due to incom- 
pleteness, their simple presence in the catalogue carries interest- 
ing implications with respect to the possible presence of this lu- 
minosity cut-off. We first consider the case that these sources rep- 
resent a smooth continuation of the universal HMXB luminosity 
function observed below ~ 10 40 erg s _1 . This luminosity is the 
Eddington luminosity (Le) for a 100 M© black hole, or 10 Lb 
for a 10 Mq black hole. One way in which higher luminosities 
may be explained is to invoke larger black holes. If there is no 
upper limit to the HMXB X-ray luminosity, then their luminos- 
ity function no longer provides evidence for an upper limit to the 
HMXB black hole mass. The situation is complicated by the possi- 
bility that some XRB s appear to radia te at super-Eddington rates 
(e.g. GRS 1915+105; iDone et alj |2004) . so an upper mass limit 
could still be pres ent without manifesti n g itse lf clearly in the lu- 
minosity function. Zampieri & Roberts] ( 2009h independently ar- 
gue that even in low metallicity environments it is not possible 
to produce black holes significantly more massive than ~ 90 Mq 
through standard stellar evolution. However, a smooth luminosity 
function would nevertheless remove an important observational re- 
sult supporting such an upper mass limit. 

The lack of a high luminosity cut-off would also raise interest- 
ing questions on the nature of the relationship between the total X- 
ray luminosit y from HMXBs in a gal axy a nd its star form a tion rat e 
presented by Gilfanov et al.l J2004ah and iGilfanov et all |2004b). 
They argue that an upper limit to the HMXB luminosity function is 
required to produce the linear regime of their relation at high star 
formation rates, based on the non-trivial statistics of number dis- 
tributions that follow a power-law. This linear regime requires that 
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the probability distribution of the total HMXB X-ray luminosity 
is Gaussian, which may only occur for galaxies with high enough 
star formation rates such that there are a fair number of sources at 
or near the upper limit. These sources dominate the total HMXB 
X-ray luminosity and remove the discrepancy between the mean 
and the mode values for this quantity which is otherwise present 
for populations drawn from luminosity functions with exponents 
1 < a < 2, and causes the non-linear regime at low star forma- 
tion rates. Without an upper luminosity limit, this non-linear regime 
should continue to higher star formati on rates. However, the data 
do not appear to continue this trend dGilfanov et al.ll2004bl) . which 
would be an important discrepancy to address. We note briefly that 
the data presented in those works do not account for intrinsic ab- 
sorption and obscuration, which may play an important role. Un- 
fortunately, attempting to account for this invariably leads to model 
based assumptions and uncertainties. 

We now also consider the case that there is a cut-off in the 
HMXB luminosity function at ~ 2 x 10 40 erg s _1 . As argued pre- 
viously, although it may be that the luminosities of some of the 
sources included in the catalogue that appear in excess of this limit 
have been poorly estimated, and some sources may turn out to be 
background quasars, this is unlikely to be the case for all 43. We 
would then be observing a population of sources radiating above 
the maximum HMXB luminosity. These must represent a differ- 
ent class of black hole, with a different formation mechanism. 
There is speculation that galaxies in the highest star formation 
regime may produce a population of IMB Hs through black hole 
mergers in dense stellar cl usters (see e.g. Gilf anov et al J 1 2004b ; 
Porte gies Zwart et al.ll2004l) . which could be seen at luminosities 
in excess of 10 40 erg s _1 . Fig. 10 in that work shows the effect 
such a population might have on the universal HMXB luminosity 
function; the cut-off instead becomes a step-like feature due to the 
upper HMXB mass limit, and above this step a far less numerous 
population of IMBHs are observed. Given the links between ULXs 
and star formation, it might not be unreasonable to argue that by fo- 
cusing on galaxies that host ULXs, we are selecting galaxies with 
high star formation rates. The presence of such high luminosity 
sources in our catalogue therefore might not be surprising even in 
the presence of an upper limit to the HMXB black hole mass, and 
may be evidence for a rare population of IMBHs. These intriguing 
sources are investigated further in Sutton et al. (in prep.). 

3.2 Specific ULX Frequency 

Here we again make us e of the complete sub-sample to mimic 
the analysis presented bv lSwartz et al.l ( I2008h . and investigate how 
the specific ULX frequency 5™, i.e. the number of ULXs per unit 
galaxy mass, evolves as a function of galaxy mass. This provides 
additional information on the type of galaxy, and hence the envi- 
ronments in which ULXs are preferentially formed. In order to es- 
timate galaxy mass for the galaxies with complete observations, 
we make use of the relations between optical colour an d stellar 
mass-to-light ratio, M/L, presented in iBell et al.l d2003l) . which 
should hold for both elliptical and spiral type galaxies (we are as- 
suming that stellar mass is directly proportional to gal axy mass). 
The d erivation of these relations utilises SDSS data dYork et al.l 
2000) and the PEGASE ga laxy evolution models (updated from 
those originally presented in lFioc & Rocca-V olmerange 19971), and 
assumes a scaled version of the Salp eter Initial Mass Fu nction 
(the same as adopted in earlier work bv lBell & de Jong 200 1]). We 
choose to calculate M/L in the B-band, and utilise the B-band 
magnitude and the (B — V) colours provided in RC3 to estimate 




L x (10 ergs g" ) 



Figure 8. The luminosity functions for sources in the complete subset lo- 
cated in spiral galaxies (top) and elliptical galaxies (bottom). Again, the best 
fit power-law model to each is also shown. 



the galaxy mass. The scatter around the relation between M/Lb 
and (B — V) is estimated to be ~0.1 dex, so we include this in 
the estimated uncertainty of our M/Lb values. The colours and 
magnitudes used have been corrected for Galactic and internal ex- 
tinction, and for redshift, and when converting absolute magnitudes 
to solar B-band lu minosities, we adopt an absolute solar B-band 
magnitude of 5.47 dCoxll200(]|) . 

Through this method, stellar mass estimates were obtained for 
86 and 35 of the 108 and 56 complete spiral and elliptical galaxies 
respectively. Galaxies without mass estimates, which were mainly a 
mixture of elliptical and small nearby galaxies (in particular dwarf 
spheroidal galaxies), were discarded from our analysis. For each 
galaxy with a mass estimate, we counted the number of ULX can- 
didates assumed to be associated with it. As with the luminosity 
function analysis, we attempt to minimize the effects of transient 
and highly variable sources by considering only sources detected in 
the longest observation of each galaxy. Then, grouping all the com- 
plete galaxies into equally spaced decade mass bins, we calculate 
S u for each bin by dividing the total number of ULX candidates 
observed by the total galaxy mass contained within that bin. This is 
expressed in equation |6]for the ith bin containing n» galaxies: 



1 Mi Eiii Mi 



(6) 



where Ni and Mi are the total number of ULXs and the total galaxy 
mass in bin i respectively, Nj is the number of ULXs associated 
with the jth galaxy within that bin, and Mj is the mass of that 
galaxy. Fig.|9]shows the specific ULX frequency plotted as a func- 
tion of galaxy mass, as well as the number of galaxies contributing 
to each bin, the number of ULX candidates in each bin, and the 
total galaxy mass in each bin. 

It is clear from Fig. [9] that 5™ decreases with increas- 
ing galaxy mass, qualitatively similar to the trend presented in 
ISwartzetal]d2008l) with a sample compiled from Chandra, XMM- 
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Newton and ROSAT observatio ns of catalogued galaxies within 
~8Mpc ( Karachentse v et alj2004h . Phenomenologically we fit the 
data (only for bins in which ULX candidates have been detected) 
with a power-law, i.e. S u oc and obtain an index of j3 — 

0.59 ± 0.07 (la uncertainties). With this relation, it becomes im- 
mediately apparent why no ULX candidates are detected in the 
few lowest mass bins. The prediction, by extrapolation, for S u in 
the highest mass bin without a detection is 7 x 10 -4 ULXs per 
10 6 solar masses which, considering the total stellar mass con- 
tained within that bin, equates to less than one (~0.4) expected 
ULX detection. For even lower mass bins, the estimated number of 
ULX detections continues to decrease; there is simply not enough 
galaxy mass contained in the low mass bins to expect ULXs to 
be detected. The lack of detections in these bins is therefore not 
inconsistent with the trend observed for the bins with ULX detec- 
tions. We also calculate S u for spiral and elliptical galaxies sep- 
arately, as shown in Fig. [10] and also fit these data with power- 
law models. The indices obtained are f3 spira i — 0.64 ± 0.07 and 
PeiUpticai = —0.5^0.5- As is clear from Fig. [10] such a model 
does not provide as excellent a representation of the data for ellip- 
tical galaxies as it does for spiral galaxies. However, we stress that 
due to a lack of ULX detections and galaxies with mass estimates, 
as well as the higher fractional contamination, the quality of the el- 
liptical data is rather poor. This is reflected in the large uncertainty 
on the elliptical index, which is in fact consistent with zero, i.e. the 
case where S u is constant with galaxy mass. Unsurprisingly, since 
the complete sample is dominated by sources associated with spiral 
galaxies, the results obtained for this population are very similar to 
the results obtained for the whole sub-sample. 

Before discussing the potential physical origins of these re- 
sults, it is important to consider possible biases that may be present. 
By considering only sources detected in the longest observation of 
each galaxy, we hope to have minimised the influence of transient 
sources as best we can. Background contaminants may ha ve an im- 
portant effect here, but as argued in ISwartz et al. I d2008l) . they are 
more likely to be prominent in the larger, higher mass galaxies, pri- 
marily because they cover a much greater sky area. This serves to 
artificially increase the number of ULX candidates and hence S u in 
the higher mass bins, and to flatten the observed trend, i.e. decrease 
/3. If these sources are important, the values of /3 presented here 
may only be lower limits. Other biases may arise due a combina- 
tion of our ULX selection criteria and the fairly modest resolution 
of XMM-Newton. By excluding the nuclear regions of galaxies and 
extended sources there will be ULXs in the observed galaxies that 
are not included in our catalogue. Since we have defined the nu- 
clear region with approximately a constant angular size, a greater 
physical area will be excluded for more distant galaxies. This is im- 
portant as the radial di stribution of ULX candidates peak s toward 
the centres of galaxies dSwartz et al.l2004l : lLiu et al.l2006l) . On av- 
erage, the higher mass galaxies are further away, so this too will 
have a more important effect on the higher mass bins, and the num- 
ber of ULXs in those bins may be slightly underpredicted. In terms 
of the exclusion of extended sources, of which star forming regions 
may contribute a large number, both galaxy star formation rate and 
stellar mass will determine the number of ULXs missed through 
being embedded within them. It is not immediately clear whether 
this should have any preferential effect on any region of the mass 
scale. Finally, there is also the possibility that some of the 2XMM 
detections are actually a number of unresolved point sources, but 
again it is not clear whether this should have a greater effect on any 
range of the mass scale than others. However, while these consider- 
ations may have an effect on the quantitative forms of the relations 
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Figure 9. The specific ULX frequency and the best fit power-law model 
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stellar mass and number of galaxies (panels C and D) as a function of galaxy 
stellar mass for galaxies with complete observations included in 2XMM. 



presented, given the similar trend presented by ISwartze t al. (2008) 
to that found here for spiral galaxies, we can be confident in the 
qualitative form of this result. 

It is clear that the number of ULXs per unit stellar mass de- 
creases with increasing galaxy mass for spiral galaxies. There are 
two main physical differences between low mass and high mass 
spirals that are likely to be combining to give this effect. The first 
is that the specific star formation rate, i.e. the star formation rate 
per unit mass, is higher on avera ge for lower mass spirals (see 
Fig. 24 in lBrinchmann et al]|2004l). ULX s are well correlated with 
star forming regions dSwartz et al.ll2009l) . and there are typically a 
larger number of U LXs in galaxies with high star formation rates 
torimmetal .1120031) . The sec ond is that low er mass galaxies typi- 
cally have lower metallicities dLee et al.l2 006). ULXs are often ob- 
served to be located in low metallicity regions (see e .g. Soria et al 



l2005llZampieri & Robertsll2009llMapelli et al.l2009llMapelli etal 
2010), and there are a number of theoretical considerations that 
suggest the Mbh > 10 Mq black holes some ULXs are specu- 
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Figure 10. The specific ULX frequency and the best fit power-law model calculated for spiral (left column) and elliptical galaxies (right column) separately. 
The vertical panels are as in Fig. [9] 



lated to host are easier to form in these regions, both due to the for- 
mation of lar ger progenitor stars, and reduced mass loss rates from 
stellar winds ( Madau & Reesl200lLlBromml2004IVink et al.l200ll ; 
Heger et al. 2003[). This combination of higher star formation per 
unit mass and lower metallicity provides a natural explanation for 
the higher specific ULX frequency in low mass spirals. However, 
we stress that this does not necessarily mean ULXs are commonly 
observed in low mass galaxies, for the same reason that there are 
no ULX detections in the very lowest mass bins. Although we find 
there are more ULXs per unit mass in lower mass spiral galaxies, a 
large number of these galaxies often need to be observed to detect 
ULXs by virtue of the fact they do not contain much mass. 

Although the elliptical data is consistent with S u remaining 
constant with galaxy mass, we briefly comment on the possibility 
that S u might increase with increasing elliptical galaxy mass as 
formally suggested by the exponent obtained, as such a scenario 



would be very interesting. When stellar bin ary systems undergo 
supernova, they experience a velocity kick dHarrison et al. l l 19931 : 
iFrail et al .ll994l : lLvne & Lorimeil 1994b . and over the course of the 
lifetime of a galaxy some fraction of its binary systems may be 
ejected. It might be that a larger fraction of binaries are ejected 
from lower mass galaxies than higher mass galaxies owing to their 
weaker gravitational potentials. This effect might lead to such a 
positive trend in elliptical galaxies where ULXs are expected to be 
LMXBs, hence star formation and metallicity effects are likely to 
be minimal. However, we stress again that the data are consistent 
with S u being constant with stellar mass, which would be a simpler 
scena rio to explain gi ven that LMXBs are expected to trace stellar 
mass dGilfanovll2004h . 
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4 CONCLUSIONS 

Through cross correlation with the RC3 galaxy catalogue, we have 
mined the 2XMM Serendipitous Survey to produce a catalogue of 
ULX candidates. Our catalogue contains 650 detections of 470 dis- 
crete sources, and we conservatively estimate that at most ~24 per 
cent, but more likely ~17 per cent of these sources are undesir- 
able contaminants, most likely background quasars in this case. Our 
catalogue offers a sig nificant improvem ent in the number of ULX 
candidates on that o flSwartz et al. (2004). which was previously the 
largest uniformly selected, dedicated ULX catalogue published to 
date. Obviously one of the key potential uses of this resource is to 
unearth and provide basic information on sources for which more 
detailed observation and/or analysis would be both interesting and 
insightful. One such source that we draw attention to is NGC 470 
HLX1, a detailed analysis of which will be included in forthcoming 
work by Sutton et al. (in prep). In addition to highlighting interest- 
ing sources, the catalogue may also be used to directly analyse what 
we hope to be a representative sample of the whole ULX popula- 
tion. 

To undertake statistical studies of the ULX population, we de- 
fine a 'complete' sub-sample of the population taken only from 
galaxies with observations sensitive enough that we would expect 
all the ULX candidates in those galaxies that meet our selection 
criteria to be detected. The limitations of these criteria have been 
discussed in detail in i]2,5l and are not repeated here. The lumi- 
nosity function of this sample is well modelled with a single, un- 
broken power-law of form N(>L X ) oc L x °'' M±lh11 . We do not 
find any evidence for a break or cut-off in the luminosity function 
at Lx ~ 2 x 10 4 er g s" 1 , as has previo usly been reported in 
ISwartz et all d2004l) and lGrimm etail 12003) and, in the latter case, 
interpreted as possible evidence for an upper limit of ~100 Mq to 
HMXB black hole masses. However, given the possible biases in- 
troduced by e.g. non-systematic uncertainties in galaxy distances, 
we certainly cannot rule out the intrinsic presence of any such cut- 
off. Additional Chandra observations with increased sensitivity of 
a greater number of nearby galaxies would help to ad dress this is- 
sue. S imilar to results obtained in previous work, e.g. ISwartz et al.l 
J2004l) . there is also tentative evidence that the luminosity function 
of ULXs in elliptical galaxies is steeper than that of ULXs in spiral 
galaxies, which would be consistent with the assumption that el- 
liptical galaxy ULXs are primarily low mass binary systems, while 
those in spiral galaxies are high mass systems. 

We find that the specific ULX frequency S u , i.e. the number 
of ULXs per unit galaxy mass, decreases with increasing galaxy 
mass for ULXs associated with spiral galaxies, and is well mod- 
elled with a power-law of form S u oc M~°- 64±0 - 07 . Although the 
quantitative form of this relation may be affected by biases such 
as background contamination, we are confident that th e qualitative 
trend i s robust, as it is similar to the result presented in lSwartz et al.l 
(2008). This trend may be due to a combination of the increase 
in specific star formation rate a nd the decrease in metal l icity with 
decre asing spiral galaxy mass dBrinchmann et al.1l2004t iLee et all 
2006). The specific ULX frequency for elliptical galaxies is found 
to be consistent with being constant with galaxy mass, as expected 
given that LMXBs should trace stellar mass rather than star forma- 
tion (Gilfanov 2004). 

The population analysis presented here is by no means com- 
prehensive, rather a number of examples of the work that is pos- 
sible with the offered resource. We hope that the construction of 
this catalogue will encourage future authors to continue this work 
and so we make the catalogue available online with this publica- 



tion ; other interesting possibilities may be to re- visit the connec- 
tion between ULXs and star formation, and to search for further 
similarities and/or differences between ULXs in spiral and ellipti- 
cal galaxies to confirm their HMXB/LMXB nature. The successes 
and limitations of our analysis highlight the importance of utilis- 
ing both the XMM-Newton and Chandra observatories to discover 
and study ULXs. A number of the limitations are due to the more 
modest angular resolution of XMM-Newton. The improvement in 
this area offered by Chandra is essential to confirm the point-like 
nature of ULX candidates, and perhaps resolve additional sources 
that XMM-Newton could not. However, the larger effective area of 
XMM-Newton provides detections with better photon statistics, par- 
ticularly at high energies, enabling analysis that would not be pos- 
sible with the Chandra satellite, e.g. detecting high energy spectral 
curvature dStobbart et al. 120061: iGladstone et al.l2009l: IWalton et ail 



1201 lUMiddleton et alj|201 1 ; Walton et al. in prep.). 

This work also highlights the importance of serendipity within 
astronomy. Catalogues such as 2XMM are vital resources for un- 
earthing new examples and studying populations of various as- 
trophysical objects. Based on the success of our compilation, we 
anticipate that future releases of XMM-Newton source catalogues, 
will contain many more ULXs again, providing a further significant 
contribution to improving our understanding of the characteristics 
and diversity of the ULX population. 
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APPENDIX A: EXAMPLE CATALOGUE ENTRIES 

Table Al. Some basic information for example entries in the presented catalogu e, including the PGC galaxy identifi er, the ado pted galaxy d istance, the 2XMM 
source ID, the detection lum inosity, potential matches in previous catalogues jColbert & Ptakll2002t CP02: lLiu& Bregman 2005: LB05: Swartz et alj|2004l : 
SW04; lLiu & Miiabelll2005l : LM05) and flags detailing whether the source is a 'new' ULX candidate (i.e. not present in any of the above catalogues) and 
whether the source detection is included in the complete subset. 



PGC 


Distance 


2XMJV1 


Detection 


CP02 


LB 05 


SW04 


LM05 


New? 


Complete 




(Mpc) 


Source ID 


Luminosity 
(10 erg s ) 


Name 


Name 


RecNo 


Name 




Subset? 


1370 


45.77 


2609 


2.25 ±0.77 










yes 


no 


1388 


45.0 


2677 


3.27 ±0.80 










yes 


no 


2248 


120.96 


4975 


5.12 ±2.18 










yes 


no 


2388 


58.21 


5291 


2.50 ±0.76 










yes 


no 


2789 


3.0 


7278 


1.35 ±0.01 










yes 


yes 


2789 


3.0 


7278 


1.18 ±0.02 










yes 


yes 


2789 


3.0 


7278 


0.38 ±0.01 










yes 


yes 


2789 


3.0 


7278 


0.87 ±0.05 










yes 


yes 


2789 


3.0 


7337 


1.62 ±0.02 








NGC 253 ULX2 


no 


yes 


2789 


3.0 


7337 


1.63 ±0.02 








NGC 253 ULX2 


no 


yes 


2789 


3.0 


7337 


2.26 ±0.04 








NGC 253 ULX2 


no 


yes 


2789 


3.0 


7337 


2.47 ±0.06 








NGC 253 ULX2 


no 


yes 


3453 


75.29 


9158 


14.38 ±2.65 










yes 


no 


4777 


34.12 


12122 


153.03 ±7.98 










yes 


no 


4801 


30.87 


12156 


1.83 ± 1.09 










yes 


no 


5098 


65.65 


12557 


14.16 ±4.07 










yes 


no 


5193 


27.45 


12818 


1.42 ±0.70 










yes 


no 


5193 


27.45 


12821 


0.99 ±0.63 










yes 


no 


5283 


73.84 


13038 


10.74 ± 2.40 










yes 


no 


5283 


73.84 


13038 


17.58 ±9.63 










yes 


no 


5974 


9.70 


14174 


2.73 ±0.15 








NGC 628 ULX2 


no 


yes 


5974 


9.70 


14223 


1.91 ±0.17 








NGC 628 ULX1 


no 


yes 


5974 


9.70 


14223 


1.04 ±0.15 








NGC 628 ULX 1 


no 


yes 


6983 


22.88 


15962 


1.57 ±0.32 






5 


NGC720ULX7/8 


no 


no 


6983 


22.88 


15962 


1.69 ±0.77 






5 


NGC720ULX7/8 


no 


no 


7252 


74.15 


16228 


11.10 ±3.64 










yes 


no 


7584 


66.47 


16743 


13.24 ±6.90 










yes 


no 


8726 


71.21 


18486 


1.78 ± 1.59 










yes 


no 


8974 


31.25 


20178 


17.66 ±5.73 










yes 


no 


9067 


83.01 


21018 


23.64 ±3.71 










yes 


no 


9578 


68.36 


23745 


16.04 ± 15.78 










yes 


no 


10122 


18.76 


25071 


0.78 ±0.29 










yes 


no 


10122 


18.76 


25091 


50.51 ± 2.15 


1X0 4 


NGC 1042 ULX1 




NGC 1042 ULX1 


no 


no 


10122 


18.76 


25091 


29.54 ± 2.40 


1X0 4 


NGC 1042 ULX1 




NGC 1042 ULX1 


no 


no 


10175 


19.65 


25272 


0.90 ±0.20 










yes 


yes 


10175 


19.65 


25272 


0.73 ±0.26 










yes 


no 


10266 


14.57 


25719 


1.08 ±0.07 










yes 


yes 


10266 


14.57 


25703 


0.91 ±0.12 




NGC 1068 ULX2 






no 


yes 


10266 


14.57 


25703 


0.86 ±0.13 




NGC 1068 ULX2 






no 


yes 


10314 


9.1 


25840 


1.02 ±0.13 








NGC 1058 ULX1 


no 


yes 


10959 


69.36 


26594 


34.03 ± 7.33 










yes 


no 


12209 


8.6 


28959 


901 -(-nil 

^.Oi U.J- J- 


IXO 6 


NdC 1 9Q1 TTT XI 


12 


wr.r 1 ?Q 1 TTT 




yes 


12286 


3.70 


29101 


1.08 ±0.06 










yes 


yes 


12286 


3.70 


29250 


5.84 ±0.06 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


11.52 ±0.15 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


8.89 ±0.16 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


6.61 ±0.13 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


7.14 ±0.14 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


11.03 ±0.21 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


8.88 ±0.19 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


7.98 ±0.23 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


10.71 ±0.29 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


7.71 ±0.28 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


9.10 ±0.43 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 


12286 


3.70 


29250 


3.23 ±0.06 


1X0 7 


NGC 1313 ULX1 




NGC 1313 ULX1 


no 


yes 
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